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Abstract 

A frame representation is used to derive a first order quasi-linear symmetric hyperbolic system 
for a scalar field minimally coupled to gravity. This procedure is inspired by similar evolution equa- 
tions introduced by Friedrich to study the Einstein-Euler system. The resulting evolution system 
is used to show that, for some classes of scalar field potentials, small nonlinear perturbations of an 
expanding Friedmann-Robertson- Walker background with a selfinteracting scalar field source, decay 
exponentially to zero or converge to constant values. As an application of our results we discuss the 
nonlinear stability of flat power-law solutions produced by an exponential potential. It is found that 
for a certain range of the parameter p, there is stability even in the absence of accelerated expansion. 

1 Introduction 

An important problem of classical mathematical Cosmology concerns the asymptotic stability of spatially 
homogeneous and isotropic spacetimes. Within this class of spacetimes, those having a nonlinear scalar 
field as the matter model have been extensively used to model early and late times cosmological scenarios 
— see e.g. [TJ [U [3] and references therein. Recently, Ringstrom [TTJ [T2] has proved that small perturba- 
tions of the initial data of scalar field cosmological solutions to the Einstein Field Equations (EFE) with 
accelerated expansion have maximal globally hyperbolic developments that arc future causally geodcsi- 
cally complete. In particular, in |12j the case of power-law solutions due to an exponential potential has 
been considered as an application of the more general methods developed in The nonlinear asymp- 
totic stability of power-law infiation solutions has also been discussed by Heinzle & Rendall [13] using 
Kaluza-Klein reductions and the methods of Anderson |17j — the latter, in turn, inspired by Friedrich's 
analysis of the stability of the de Sitter spacetime [3|. 

A natural way to analyse the stability of spacetimes is to ask whether small perturbations of a given 
solution to the EFE asymptotically decay to a background solution. Most of approaches to this question 
have been limited to the use of linear or higher-order truncated perturbation theory, and thus, they never 
take fully into account the nonlinearity of the EFE — see e.g. [131 dSl dS] and also [22] • This type of 
analysis has been hampered by the lack of a suitable formulation of the EFE for which the theory of 
systems of first order hyperbolic partial differential equations can be applied. 

In [5] Friedrich has introduced a frame representation of the vacuum EFE. The evolution equations 
implied by this alternative representation of the equations of General Relativity constitute a first-order 
quasi-linear symmetric hyperbolie system (FOSH). In general, these systems are of the form 

A°(i, X, u)9tu - A^ (t, X, u)9jU = B(<, x, u)u (1) 

where u — \i(x,t) is a smooth vector-valued function of dimension s with domain in E x [0,T] with S 
a spacelike 3-dimensional manifold. Moreover, A°, A^, j — 1, 2, 3, and B denote smooth s x s matrix 
valued-functions. The matrices A'' and A-' are symmetric. In addition, the matrix A*' is positive definite. 
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The operators dt and dj stand, respectively, for the partial derivatives with respect to the coordinates t 
and . 

The construction for vacuum spacetimes given in [S] has been extended in [SI [7] to the case of a perfect 
fluid using a Lagrangian description of the fluid flow — see also [211 125] . In both the vacuum and the 
perfect fluid cases the introduction of a frame formalism gives rise to extra gauge freedom. This freedom 
is associated to the evolution of the spatial frame coefficients along the flow of the time-like frame. If one 
fixes conveniently this gauge (using, for example, the Fermi gauge) one obtains a hyperbolic reduction 
for the evolution equations. As a consequence, given smooth initial data satisfying the constraints, local 
in time existence and uniqueness of a solution to the EFE can be established — see e.g. [Tj [19] and also 
[H for details. 

A natural way of performing a stability analysis is to consider a sequence of smooth initial data sets u*^ 
for the EFE satisfying the constraints equations on a Cauchy hypersurface S. The sequence is assumed 
to depend continuously on the parameter e in such a way that the limit e — renders the reference 
solution ii. In particular, one can write the full solution to the EFE as the Ansatz 

u'=u + eu^ (2) 

where u*^ is a (nonlinear) perturbation whose size is controlled by the parameter e. Using the Ansatz ^ 
in equation ([ij, and writing 

B(u + eu') =B(u) + eB(u,e), 
A^(u + eu") = A^(u) + eA^(u, e), ^ = 0, 1, 2, 3 

we are led to consider the following initial value problem for the nonlinear perturbations: 

(^KP{t, X, u) + eA^{t, X, u, e)j dtu — (^A.\t, x, u) + eA\t, x, u, e)j djii = ^B(<, x, u) + eB(t, x, u, e)j ii, 

ii(a;,0) = Uo{x). 

(3) 

A particular approach to the existence and stability of solutions to the Cauchy problem (jS]) has been 
discussed in [51 [SI [H] . In this approach, the stability of solutions follows from the existence of eigenvalues 
for the nonprincipal part of the linearised system (e = 0) having a negative real part. Stability results 
for the case where the coefficients of the linearized system are constant matrices are well known [5] . In 
the case where the system is only strongly hyperbolic, the inner product in has to be replaced by the 
so-called "H-inner product — see [9] . A procedure to analyse stability in the case of systems where B has 
vanishing eigenvalues has been given in [S] — see also [101 [H] . These methods can be easily generalized 
to systems of the type considered here where the matrices B. A^ are not constant but depend smoothly 
on time — see also [5D]. 

The approach described in the previous paragraph has been applied by Reula [5D| to the Einstein- 
perfect fluid system of [6] with a cosmological constant, to prove the exponential decay of nonlinear 
perturbations for a wide class of homentropic fluids in flat Friedmann- Roberts on- Walker (FRW) back- 
grounds. The analysis of [ID] gives, as only restriction on the equation of state, that the sound speed is 
less than one-third of the speed of light. An advantage of this approach is that it avoids the problem of 
gauge-dependence in perturbation theory. 

In the present article we pursue a similar approach to the one used by Reula in |20| to analyse the 
nonlinear stability of FRW spacetimes with a nonlinear scalar fleld. To this end, we flrst construct a 
flrst order symmetric hyperbolic system for the EFE with a scalar fleld as the matter source. This 
construction is performed by splitting the wave equation for the scalar fleld into two first order equations. 
In our analysis, the scalar field is used to construct an adapted orthogonal frame. Written in terms of 
this adapted frame, the energy-momentum tensor is diagonal, independently of further gauge choices. A 
similar construction has been considered in the analysis of linear perturbations in [151 118] . 

An important difference between the system treated in the present paper and the one in reference 
|20j is that in our case the matrix B has eigenvalues which are zero. These vanishing eigenvalues appear 
in the block containing the perturbed matter variables and expansion. As it will be discussed in the 
sequel, this feature is related to the fact that a subspace of the perturbations converge asymptotically to 
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constant values. Our analysis of the eigenvalues of the nonprincipal part of the linearized system uses 
results from the so-called Routh-Hurwitz problem. In particular, for an ever expanding FRW-nonlinear 
scalar field background with spatial topology T"^, i.e. ([0, oo) x T'^,g^(j)), we show that, for scalar field 
potentials satisfying certain inequalities, nonlinear perturbations exist and asymptotically exponentially 
decay to zero or converge to constant values. In the suitable limit, we recover the results of Ringstrom in 
[llj . As an application of our results, we show that for power-law solutions produced by an exponential 
potential with 

P > ~ 0-454 

and flat spatial sections, sufhciently small nonlinear perturbations decay exponentially to zero or con- 
verge to constant values. We have confirmed this behaviour in some simple numerical computations of 
linear perturbations. Our analysis shows that nonlinear stability of power-law solutions with exponential 
potentials is possible even in the absence of accelerated expansion — a possibility that, hitherto, has not 
been discussed before. 

Our main result is proved in the Sobolev norm H'' (T"^) for fc > 5. As a consequence of the use of 
results from the Routh-Hurwitz problem, the conditions we obtain on the potential are, in a sense, sharp 
and cannot be improved by the present methods. Given a particular background Cosmological solution 
to the Einstein-nonlinear scalar field system — calculated, say, numerically — one can verify explicitly 
whether the conditions are satisfied or not. 

Overview of the article 

The paper is organized as follows: in Section 2, we recall Friedrich's frame formulation of the EFE. In 
Section 3, we discuss some relevant properties of scalar field satisfying a nonlinear wave equation. In 
Section 4 we discuss the conditions under which the Einstein-Friedrich-nonlinear scalar field system is 
well-posed, in the sense that it forms a symmetric hyperbolic system — see Thcorcm[TJ Finally, in Section 
5 we give the conditions for which there is convergence to constant values or exponential decay of small 
nonlinear perturbations on a FRW-nonlinear scalar field background. This is the main result of the paper 
and we summarise it in Theorem [51 We use units such that SttG = c = 1. 

2 Friedrich's frame formulation of Einstein Field Equations 

In this section we provide a brief introduction to Friedrich's frame formulation of the Einstein field 
equations. The basic equation of Friedrich's construction is the contracted Bianchi identity. From the 
latter it is possible to deduce hyperbolic propagation equations for the conformal Weyl tensor for a wide 
spectrum of gauge choices. 

2.1 Basic definitions 

Let (A^,g) denote a 4-dimensional spacetime with g a Lorentzian metric of signature (—,+,+,-!-). In 
order to implement the frame formulation of the Einstein field equations, one defines locally an orthonor- 
mal moving frame or tetrad with respect to the metric g in an open neighbourhood lA <Z Ai. The frame 
is a set {ea] of linearly independent vector fields in the tangent space Tp (A^) at each point p E U such 
that 

g {ea,eb) ^ r]ab, a, fe = 0,...3, (4) 

where rjab ~ diag (—1, 1, 1, 1). The norm of a vector field, v E Tp {A4), in an orthonormal frame is defined 
as 

\vf = g{v,v) ^ v'^v'Tjab 
and in terms of a coordinate basis set {da} we have Ca = e^d^. Condition (g]) gives 

Vab — C-a 9 fiu ■ 

where = 0, 1, 2, 3 denote spacetime tensor indices. The frame commutator is written as 

[ea,e6] c^^bSc, (5) 



3 



where c'^^jj are the structure coefficAents. The dual basis or coframe is the set of hnear forms {9'^} belonging 
to the dual space T* {^A) at each point p defined by the pairing {6^ ^ e^) = 5\. In terms of the dual 
basis we can write condition (|4]) as 



9 = -(^")'+E(n^ 



The spacetime (Levi-Civita) connection in an ortho normal basis is defined by 

where 7^^,^^ are the connections coefficients. The covariant derivative of a tensor in M can be written as 

^1 fa"gi---gs 1 qia"f---qs 1 q,a"gi---f 

The torsion free and metric compatibility conditions imply, respectively, that 

The equations for the frame coefficients {e,/'} are given by equation ([5]) in terms of the connection 
coefRcients. In turn, equations for the connection coefficients are obtained from the Ricci identity 

R%cd = ec - e, i^bc) + ^fclL - iWbc - 7V (7^ - ■ (6) 

The Ricmann tensor can be decomposed in terms of the conformal Wcyl tensor C and the Schoutcn 
tensor S as 

R°'bcd = C^bcd + <5°[c'S'd]6 — VblcSd]"^- (7) 

For future use we introduce the Friedrich tensor F via 

Fabcd = Cabcd " Va[cSd]b: (8) 

and its dual with respect to the last pair of indices 

* Fabcd = *Cabcd + -^Spbe^acd, (9) 

where eabcd is the usual Levi-Civita totally antisymmetric symbol with £0123 — 1- In terms of the Friedrich 
tensor, one finds that the contracted Bianchi identities read 

Vancd=0, Vanned = 0. (10) 

2.2 Orthonormal decomposition of the field equations 

The field equations of Friedrich's frame formulation of the Einstein equations arc given by equations ([5]), 
dl]) and pH]) . together with the decomposition ([7]). The independent variables of the system are 

(Sa'^, 7"bc, C'^bcd, Sbc) ■ 

In what follows, we will shall decompose the field equations and relevant tensors in terms of their 
parallel and orthogonal components with respect to the time-like frame. We write TV = eo and set 

AT = N'^ea, = Jo", 

where Na ~ —Sj' in our signature. In terms of these objects, tensor fields which are orthogonal to the 
timelikc frame- vector are defined by 

Ta,...a,...a,N''' =0, p = 1, 2, . . . , 

Next, one defines 

hab = Vab + NaNb, 
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where the projector onto the orthogonal 3-subspaces satisfies ha'^ = Tf"^hab- The spatial covariant deriva- 
tive is then given by 

n T — h ■ ■ - h P^\Ii,T 

In particular, one has that 

Dahd = 0, Daebcd = 0, 
where ebcd is the spatial Levi-Civita symbol and the indices run from 1 to 3. 

In order to further proceed with the geometric decomposition one defines the acceleration vector by 

a = Voeo = 7^ooep, p= 1,2,3. 

It follows then that qP = j^qq or equivalently, Up = I'^po- We will also consider the so-called Weingarten 
map given by 

X{ea) ^VaCo = YoaSp, a,p= 1,2,3. 

It follows then that ~ J^oa- The tensor Xa6 can be written in terms of its irreducible parts as 

Xab = 7°6a = {X^'^)ab + ^X^ab + {x'^)ab, 

where {x^'^)ab, X^ {x"^)ab denote, respectively, its symmetric trace-free, trace and antisymmetric parts. 
If the flow of Bq is hypersurface orthogonal then one has that {x^)ab ~ and that 

^£Nhab ^ X{ab} ^ {x'^'^)ab+^Xhab, (11) 

where £n denotes the Lie derivative along TV and WaN^ = —Nao'^ + Xa' ■ 
Finally, the 4-dimensional Levi-Civita symbol is also decomposed using 

e-abcd = 2ea6[c-^d] " ^^^la^^cd- 

Now, defining Fbcd = aF°'bcd it follows that the first contracted Bianchi identity can be written as 



Fbcd = Nb 



F^Nd - F,5do^c + 2F,'o[,iVrf] - NbP;,,, + Fl,, = 0, (12) 



where contractions with N are denoted by the index 0, and the prime ' indicates that the remaining 
indices are spatial. For example, F'bod = hb'^N^hd'^Fqrs- 

Given the vector TV, the Weyl tensor is uniquely determined through its electric and magnetic parts. 
These are defined, respectively, by 

Eab = haPhbrn-N'^Cpgcd, Bbd = hb^hd'^N'^N-Clp,^. 

In terms of the latter, the Weyl tensor and its dual can be written as 

Cabcd = 2 [la[cEd]b - h[cEd]a\ - 2 [TV[ci3(J]pe''ah + N[aBb\p^^ cd\ (13) 
*Cabcd = '^N[aEb]p(-^ cd — 4i?p[ae6]'' [c-^d] — ^^[aBb][cNd] ~ Bpqt'^ ab^'' cd (14) 

where lab = Kb + NaNb. 

3 Nonlinear Scalar fields in the frame formalism 

In this section, we introduce a description on nonlinear scalar fields which is particularly well adapted to 
the analysis to be carried out in the present article. 
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3.1 Basic equations 

In what follows we will consider scalar fields (f> £ C°° {M.) with a smooth positive selfinteracting potential 
V ((/<), and energy-momentum tensors of the form 

r = tA®V- Ql^l' + V(0)^g, 

where we have defined the 1-form 

with V the spacetime connection. Accordingly, we define 

V'a = V' (ea) = (15) 

where we have written 

■0 ^ -00 = £N't) (16) 

and 

^'^ = hj'i^b = Da(l). (17) 

The components of the energy- momentum tensor T with respect to the tetrad {eq} are given by 

Tab^^'a^b- Q|t/'P+V(0))77a6, (18) 

while its trace is 

r = -|V'p-4V(0). 

The Einstein field equations in the form 

R^T \g (TrgT), 
then imply for the components of the Ricci tensor that 

Rab = V'aV'6 + V(0)77a(,, 

while the Ricci scalar is given by 

R=-T= +4V(0). 

From these expressions it follows that the components of the Schouten tensor with respect to the frame 
{ca} are given by 

3.2 Gauge considerations 

In order to construct and adapted frame to our particular problem, we let = aeg. It follows that 

r = a6^, (19) 

so that 

a = --0 and D^'if) ^ 0. 

Accordingly, 



|0P =g(-0,0) = «2^oo --a', « = ±v/-|0P. (20) 

If the vector tp is taken to be future oriented, then one must choose a to be positive. In terms of a 
coordinate basis the latter implies 

0M = aeo^ = _0eo^ gq'' = (21) 

and 

Da(j) - 0, e;^V,,0 = 0. (22) 
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With this choice, we have that 
and therefore 



Tab = + V(0) j NaM + Q^' - V((/.) ) hab, (23) 

1/5,2 A - 1/1 



Sab - 3 (^2^ - "^'('^) j + 3 [2^ + "^('^) ) ''-b. (24) 

Using (|15p and (j20[) , the equation for the conservation of the energy- momentum tensor takes the form: 

y^'Tab = (^i'^NaNb + QV^2 _ (25) 

= 2^NbN'' (VaV) + V^' {Nb i^aN'') + TV" (VaA^b)) + Vb f _ (26) 



2^j£^i, + + ij^) Nb + tP^ab + i>Vbi' - 0. (27) 
d(j) 



From the latter, projecting with respect to the timehke frame one obtains: 

iv' {V'Tab) = 0, + + ^ = 0, (28) 

{V'Tab) = 0, D,^ + = 0. (29) 

Moreover, using that Da<i> = in the orthogonal subspaces to Bq, one obtains from equation ([Sj that 

K,e;,]0 = 2(x^)afc^ = O, 

which implies that 

{X^)ab = 0. (30) 

Remark. Following Friedrich in [7], one could as well have defined 

VT^f, = gfc + qNb, Jab = V[„gfc]. (31) 

Then, instead of using the condition on the vanishing of the divergence of the energy-momentum tensor, 
one could include the equations q = Q and = as a part of the equations determining the Einstein-scalar 
field system in the frame representation. Once the gauge is fixed, the first equation in (|3ip appears in the 
reduced system of evolution equations while the second part is regarded as a zero quantity — see equation 
(4.44) in [7|. It can be shown that the zero quantities satisfy a system of subsidiary evolution equations. 
For this, it can be shown that the zero quantities vanish if they are zero on the initial hypersurface. For 
the quantity qb the relevant subsidiary equation is given in equation (4.70) of [7j- We also notice that the 
evolution for the acceleration can be computed from the tensor Jab- 



4 The Einstein-Friedrich-nonlinear scalar field system 

In this section, we derive a first order symmetric hyperbolic system for the EFE coupled to a nonlinear 
scalar field. Making use of the Bianchi identity and the energy- momentum tensor of equation p8|l . we 
derive the propagation equations for the electric and magnetic parts of the conformal Weyl tensor. After 
fixing the gauge, we complete the reduced system of evolution equations by deriving equations for the 
frame and the connection coefficients. In the last part of this section we make some remarks concerning 
the hyperbolicity of the system. 
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4.1 Basic expressions 

Wc start by computing the various components for the Friedrich tensor F. Using equations ([T^ and 
one finds that the definition Q imphes 

Kbcd = Bbpe^cd = -Kbdc = ^Hocd: 

Fobod = Em + 1 Qv-' + V(0)^ /^brf = -J^^,,o, 

^I&cd = "2 {hb[cEd]a - ha[cEd]b) - ^ QV"^ + V(<?!))^ (/lacft-db ^ /larf/lcb) , 



with the nonvanishing traces 



h-^h'"'F:,,,d = -^^'-m- 



(33) 



Using the expression © with equations (|14p and (HH), we get the following components of the dual *F: 

* Tp' ^ Tp' n + 17^/ ^ jpf n 

* Tp' * rpf D 

'^^locd = ^ape^cd - ^ - Vi^)^ eacd, ^^^^ 

''FUcd = -EbpePcd - I + V(0)) efc,rf, 

* it' 



FobOd — ^fcdi *Fabcd — BpqS^ ab^^ cd- 

4.2 The Bianchi equations 

If one substitutes the expressions for the Friedrich tensor derived in the previous section into the first 
Bianchi identities ()12p one obtains the following relations for the components of the zero quantity Fabc'- 

FqcO = ^ "^-iv^oOcO + ^'-^qOcO + Xc Fqq^q — XFqqcO — X"^ (^(}bcO + -^qOcb) + (^ObcO + ^OOcb + ^bOco) ' 
-^Ocd ^ ^ "^AT-Fooc^ + D''F^Q^d + ^''-^Obcd + ^'-^qOcd ~ X'^'^ F'^f^^d ~ XFoocd ~ X'^c-^gOOd ~ X'^d^ijOcO 

+ Xc ^OOsd + Xd -^OOcs + ^c^oOOd + '^'i-^OOcO' 
Fbod = ^ "^AT^ohod + E'^F'^f^Qd ^ XFobod ^ X^b^aOOd ~ X°'^F^f^^d + Xfc^'-PosOd + ^b^oood + Xd^^ObOs 

+ ^"^F^ijQd + a'^-Pobcd' 

^bcd ^ "^Jv^Obcd + E"" F^i^^d + ^''^gbcd + (^b^oocd + ^cFgi^Qd + Od-Fobco) ~ X-^Obcd 
Xb -^Ogcd ^ Xc -^Obqd ^ Xd ^Obcg X b^ qOcd X qbOd X d^ qbdCn 

(35) 

where we have used the fact that is anti-symmetric in the last two indices — see e.g. [S]. Similar 
relations hold for the dual 

Remark. In reference [7] — cfr. equation (4.47) — suitable zero quantities are defined by using the 
decomposition in terms of irreducible components of F . 
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4.2.1 The evolution equation for the electric part of the Weyl tensor 



An evolution equation for the electric part of the Weyl tensor can be obtained using the third equation 
of ([551) together with the expressions ([5^ and and then symmetrising with respect to the indices 
(bd). One obtains the equation 



F!, 



1 /I 



{b\0\d) 



pa 



+ 2aaBp(^bedr - ^xEm + 2x\bEd)q + Sxib^E^), - hdbX""Eac - 3 (^'' - V(0)) X{bd)- 
Similarly, using equation we get 

-^(b|o|rf) = ^ £NEbd + DaBp(^d'^hf'^ + '2aaBp(^hedf^ — 2xEm + 2x'^{bEd)q + SXib'^E^-^g — hdbX'"^Eac 



The trace of the previous expression is given by 

h'^'Fl,\o\s)--li^'x-l£. 



^2 + v(0) 



which is the evolution equation for the scalar field — i.e. the equation expressing the conservation of 
energy. The evolution equation for the components of the tensor Eat is obtained by taking the difference 
of the last two equations so that: 



£ mEl 



N-t^bd 



DnB. 



a^p{d^b) 



-~2a„B 



1 



2xEm + 2x\bEd)g + iXib'Ed^g - hdbX""Ea 



7;^ X(bd) - i^xh 



bd 



Finally using (j30p we can write the last equation as 

£.Ebd - DaBp^debf'' = (x'^^)m - \xEbd + 5 (x^^)f, Ed)^ + 2aaBpi^bW 



(36) 



4.2.2 The evolution equation for the magnetic part of the Weyl tensor 

An evolution equation for the magnetic part of the Weyl tensor can be derived from the third equation 
of ((55)) using the expressions ([Ml) . A computation yields 

FbM = - £mBm + [2Ep^beafd + \ QV-' + V(0)^ e^a^ - xBbd + X^bBad + 2X(6"Bd)a 

+ 2a,Bpbe''Pd + a^Ep^ePM + x'"" Bp^e^ ^b^'' cd- 

Now, since Bm is a symmetric tensor, all the information about its evolution is contained in the sym- 
metrised expression F^b\Q\d)- Consequently, symmetrising the previous equation with respect to the spatial 
indices {bd) we find 

£NBbd - DaEp^bf^d)"^ = -'2aaEp(b<^df"' + '2X{b"'Bd)a + X°'{dBb)a - XBbd + XacBpqt^"' (b'^df" , 

so that using ([50)) we have 

dtBbd - DaEp^b^df^' = -^-xBbd + 3 [x^^y (bBd)q + (x^^),, Bp^eP^b^d)'' ~ 2aaEp(bedf\ (37) 
The principal part of the evolution equations p6p and ([37p is given by the symmetric matrix 



/ eo 











--2D1 


hD2\ 





eo 










-5^3 








eo 




5^3 












eo 


















eo 







-5^3 











eo / 
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4.3 The Lagrangian description and Fermi transport 

In order to deduce the remaining evolution equations, we will adopt a Lagrangian description. This point 
of view amounts to requiring the timelike vector of the orthonormal frame to follow the matter flow lines. 
Accordingly, we introduce coordinates (t, x) such that 

eo^dt, eo'' = S^. (38) 

This particular choice is equivalent to setting O'^ ~ 9''jdx^ while at the same time fixing the lapse function 
to onc0. With this choice (since £n = dt), we have from equations (jl6p . (|28|) . and ([^^ that 

dtct) = ^J = N^i^a = -a < 0, (39) 

at^ = -V'x-^, (40) 

elV = -e'JD,^. (41) 

Now, the timelike coframe is given in terms of the natural cobasis through the relation 

0" = dt + /?jda;^ 
while the spatial frame vectors are found to be 

e', = {ea')-'{d,-P,d,), e',° = {ejr'/3„ e'J = {9a=r\ (42) 
It then follows from equations (|4T|) and (|42)) that 

ft = -^D,cp. (43) 

Thus, since Pj is nonzero, the surfaces of constant time are not necessarily spacelike for the characteristic 

cone and this could be a problem for the hyperbolicity of the system — see [MJ [35] . 

Finally, the remaining frame components are chosen to be Fermi propagated along eg. That is, we require 

VoGq - {g (ca, Voeo) eo - g (e^, cq) VoGo) = 0. 

The latter condition implies 

7'%o - 0. 

4.4 Evolution equation for the frame coefficients 

As already mentioned, the evolution equations are obtained from the relation ([5]) to yield 

[eo, e^] = flbeo - T'^obe^, 

where it has been used that 7'^6o = (Fermi gauge). Therefore, the evolution equations for the remaining 
frame coefficients read 

a / z c I i 

"tCt = -Xb , 
dte',' = flb - Xb%'. 

which, together with (|42D . imply propagation equations for the components of the metric in the local 
coordinate system. In particular, one has that 

with Pj given by equation (031) — see also equation (6.2) in [53] for arbitrary lapse U. Finally, since from 
(|30| we have {x^)ab = 0, the evolution equations for the frame coefficients can be written as 



(44) 



^Sce also 1241 . where a symmetric hyperbolic system was obtained for the Einstein-Euler system. This construction holds 
for an arbitrary Eulcrian picture. 
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4.5 Evolution equations for the connection coefficients 

The equations for the connection coefScients arc obtained from the splitting of the Riemann tensor with 
respect to the frame {ea}. In general, we have that 

Fd^bod = eo {xdb) - Ddttb - ttbttd + XpbXd^ ~ XdpYba ~ XpbYdo, , , 

R'%cd = D,Xd^-DdXc''-a^{Xcd~Xdc), ^ ' 

n bed — n bed + Xc Xdb — Xd Xcb — 1 bO [Xcd — Xdc) , 



where R°'bcd denotes the Riemann tensor constructed only with the spatial connection coefficients y^ab- 
The first two identities give evolution equations once the Lagrangian gauge is introduced. The remaining 
two equations are the quasi- constraints for the connection coefficients — see [211 No equations for the 
connection coefficient associated to the acceleration can be deduced from these identities. In the sequel, 
it will be shown how evolution equations for the acceleration can be obtained for our particular problem. 

From equations (|45p. we can also deduce two important equations relating the Ricci tensor to the 
connection; 

Rqq = -eo (x) + Dpa^ - Xp''Xb^ + apa^ 
R'od = DcXd'~DdX-2a-{x^)^^. 



The first identity in P5|) . together with the conditions for the Lagrangian and Fermi gauge provide 
the equation 

da\d = -7^pX/ + 2h''''Xd\j,ab] + BdpeP\, 

describing the evolution of the spatial connection coefficients ^"^ab- To obtain the last equation it has 
been used that R"^bM = C""bod = Bdpt'^°'b- Moreover, since from (pO| we have {x^)ab = 0, we can write 

dt^bd = -\xi\d + \xh''Phd^ab} + BdpeP\ - {x'^)/l\p + 2/.°^ {x'^) ,^^aby (47) 

The evolution equation for the part of the connection described by Xbd is obtained from the second 
identity in In order to do so, first, we will derive the evolution and the quasi-constraint equations 

for the acceleration. The evolution equation for the acceleration can be obtained from 

[eo, e^ V = c^oceo{i>) + cPoce'p(i/') 

= 7%(at^) + (7^o-7V)(i?P^) 
= a, {dti') - Xc" {Dp^) , 

where the Lagrangian and Fermi gauge have been employed. Now, expanding the left hand side and 
making use of the evolution and the quasi-constraint equation for the energy-momentum tensor of the 
scalar field one has that 

dtac - DcX = -Xc a-p + [x + 



tp d(j) 



so that using the second equation in (|46p we arrive at 



dtac - DpXc" = -2 (x^)^, aP - x/a^ + ( ^ + ) a,. (48) 



d(j) 



Finally, since for our particular problem one has that (x )ab = 0, this last equation takes the form 

0,ac - Dpxc' -(-^ + Ix) ac - (x^^)/a,. (49) 



d(j) 3 

In the case of the quasi-constraint a computation yields 



Dcab~Dbac^2[x+^^]{x'')cb- (50) 
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Thus, making use of this equation in the second identity in equation (j45p and recaUing the properties of 
the Fermi gauge one finds that 



1 



1 dV 



ifi d(j) 



Idb 



where it has been used 



R'°bod = C'\m + t^S'^^Sm + T^hbdSo" = Em + - (V(0) - V') hM. 



The principal part of the combined system of equations (|48p and ([5T|) is given by 



/ eo 




-D2 


-^3\ 


-A 


eo 













eo 













eo / 



3 



(51) 



Since {x^)ab = 0, then Xbd is symmetric and in this case, we have after symmetrisation that 

2dtX(bd)-'^D(,a,) = ^ (^V(</)) - ^2 _ 1^2^ /ifcrf-lx(x^'^)b<i-2(x^'^)/(x^'^)p, + 2abad + 2£;bd. (52) 
Also, from the first equation in (|46p it follows that 



9tX - = V(0) - V' - 



,ST\ 



where (x'^^)^ = (x'^^)a6(x'^^)''^ and = Uca". 
4.6 Hyperbolicity considerations 

The system consisting of equations dM]), dSZ]), dH]), (HZ]), dMl) and ([51]) can be written matri- 

cially as 

A°dtu - Af e;(u) = B(u)u. (53) 

As discussed in [21], these systems are not hyperbolic in the usual sense as, in general, the time lines are 
not hypersurface orthogonal and the "spatial" frame vectors have components in the time direction 
— see (|35|). Since the surfaces of constant time t are not necessarily spacelike, this type of system is 
referred to as a quasi FOSH system — see [53]. In terms of the partial derivatives, equation ([5^ reads 



with 



A"(u)afU - A^{u)djU = B(u)u 
A°(u) = A°- APe;^ A^ (u) = APe;^'. 



(54) 
(55) 



In order to have a well posed initial value problem, the matrix A*'(u) must be positive definite. This 
is the case as long the quadratic form 



&=1,2,3 



Aft- 



(56) 



is positive definite - see Proposition 9 in |24) . In the next section we will consider a reference solution 
admitting a foliation by homogeneous spacelike hypersurfaces. As a consequence, the linearisation of the 
system ^ is well posed without the need to control the smallness of ft. The smallness of these terms is 
taken care by the perturbation fields — see also [50]. 
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Written in terms of partial derivatives, our system of evolution equations reads 

dt(j) = tp, 



dV 



2d, 



'tX{bd) 



2e' 



dtEi 



bd ■ 



<^''"'{b\e'a"dtBp\d) 



dtBbd - e'^P 



(de\adt\Eb)p 



e'^\de\Jd,\Eb)p 



dte'b' 
dtl'\d 



v(0) - r 

ST\ 



1 



j hbd - il^bd + Ydb) CLp 
jx ix'^bd 2 (X"^)/ (X"^),, + 2abad + 2Ebd, 

ST\ 



f 2dV 2 



ab 



l'c,{x'^)% 



7 



{x'^),d-\xEbd + h{x'%Ed), 



+ 2aaBp(^i,edf - l'^ paB q(d<^bf 



<^^°'{dl'^b)aBpq, 



--^xBbd 



iix'^MB 



(b ^d)q 



'2aaEp(^b<^df°' - I'^paEq^b^d) 



~ ^^'\bl'd)aEpq, 



;Xeb 



[xnb%\ 



1 



-^Xe',"-(X^^),V.° 



ab, 



3 

-lxi'\d + Ixh'^^hdipab] + BdpeP\ - (x^^)/7'" 
+ 2h'^P{x^^) d[pab]. 



bp 



This system has clearly the form given by equation f[M|) . If one writes 

= V, z^, w^, x^, y^) , 

where 

= (Xll, X22, X33, (X^^)i2' (X^^)i3: (X'^^)23' "l' ^2, 03) 
W"^ = (£^12, Ei3, E23, B12, Bi3, B23) , 

= (ei°, 62°, 63°, ei\ 6l^ 6l^ €2^ 62^ e2^ 63\ 63^, 63^) 

= (7^22, 7^33, 7^23, 7^iii 7^33, 7^31, 7^11: 7^22, 7^12) , 
then the matrices given in equations (j54p and (j55p have the explicit form 



A"(u) = 



with 



/I2x2 
















02x2 
















V 


■^9x9 






AO 

■^6x6 







I21X21/ 




AJ(u) = 





^ 


A^ 

-^9x9 







A^ n 

■^6x6 ^ 

O21X21 






/ 1 

















p 
















1 

















6 

^2 

















1 

















-6 
'^3 

















2 








6 - 
^^2 


6,0 







-^9x9 — 




















2 




2 


6 






6 

'^3 


-6,° 

-e2° 
















-6 

62 


-6 





1 

















-6 

^2 





-6 





6 

63 





1 











V 





p 

63 





-6 _ 

^1 


6 

^2 








1 J 
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A" - 

^6x6 ^ 



/ 1 











_ip 
2^1 


ip \ 
2*^2 \ 





1 





ip 





-ip " 
2 '-3 








1 


_ip 
2*^2 


ip 
2'^3 








ip 


_ip 
2^2 


1 








2^1 





ip 





1 





I ip 

\ 2^2 


-ip " 











1 / 



-^9x9 



and 



A^ 

-^6x6 





/o 



















o\ 
























62' 






























p ^ 

''3 





















^2 


e,^ 


























p ^ 

^3 





























p ^ 

''3 













f^2 


p ^ 





















62' 


p ^ 





p ^ 

^3 















V 


e,^ 





p ^ 


p ^ 

^2 














/ 


















i 


-ip A 

2'^2 

















-ip ^ 

2*^1 







ip ^ 
2*^3 

















ip ^ 

2*^2 


-ie 


j 

3 












-ip 

2^1 


ip 

2^2 



















ip 





-ip ^ 

2*^3 















V 


-ip 

2*^2 


ip 

2'^3 


















/ 



It can, therefore, be verified that A''(u), A-' (u) are symmetric and that, furthermore, A''(u) is positive 
definite as long as ([55]) is satisfied. We summarise the results of this section in the following theorem: 

Theorem 1. The Einstein- Friedrich-nonlinear scalar field (EFsf) system consisting of the equations 
in (|57p forms a quasi-linear first-order symmetric hyperbolic (FOSH) system for the scalar field, its 
momentum- density, the frame coefficients, the connection coefficients and the electric and magnetic parts 
of the Weyl tensor, relatively to the slices of constant time t, as long as the quadratic form 



E 



=1,2,3 



-0 1/) 



is positive definite. 

Using the standard theory of symmetric hyperbolic systems one can then conclude the local existence 
in time and uniqueness of smooth solutions for the evolution equations implied by the Einstein-nonlinear 
scalar field. In order to conclude the existence of solutions to the full Einstein-scalar field equations one 
has to verify that the constraint equations are satisfied during the evolution if these hold initially. This 
will not be discussed here as the purpose of the present work is to study the evolution of the perturbations. 
In any case, it follows by general arguments that these are preserved during evolution — sec e.g. pOlfTOl IT]. 

Remark. As part of the procedure to make the system ([57| explicitly symmetric hyperbolic one has to 
divide the evolution equation for the acceleration (05]) by ■ip'^. This could imply that the system is not well 
behaved when -0 = 0. An inspection shows that the potentially troublesome term is the one containing 
the first derivative of the potential, which, by virtue of the EFE must be zero in the asymptotic limit. 
Thus, the coefficient 

V 

(60) 



lim — 

t— >-oo if) 



const. 



must be required to exist. Moreover, we notice that our gauge restricts the set of solutions, since it 
only admits scalar fields </> which are monotonically decaying functions. Accordingly, the stability of 
background solutions for which — >■ cx3 cannot be studied by our method. 
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5 Stability Analysis 



In this section we use the symmetric hyperbolic system derived in last section to show that, for some 
classes of potentials, the evolution of sufficiently small nonlinear perturbations of a FRW-nonlinear scalar 
field background prescribed on a Cauchy hypersurface with the topology of a 3-torus T'^, either converge 
to constant values or have an asymptotic exponential decay. 



5.1 The background solution 

As it is well known, the metric of a Friedman-Robertson- Walker (FRW) spacetimc — i.e. an homogeneous 
and isotropic spacetime — can be written as 

ds^ = _di2+ j^^y J.^dzMx^ 

where 

and the constant fc = — 1, 0, 1 is the curvature of the spatial hypersurfaces. Since the metric is conformal 
flat, it follows that 

Ebd = Bbd ~ 0. 

Furthermore, if one foliates the spacetimc with the surfaces of constant t, one has that 
For such background metrics 

X = 3a/a = 3H, 

where ' denotes differentiation with respect to time t, and H is the so-called Hubble function. In the case 
of a FRW cosmology the Einstein-scalar field system reduces to the evolution equations 

subject to the Friedmann-scalar field constraint equation 

H^t)-U^{t)-\vi^) = -\. (62) 
6 3 

Now, the gauge conditions for the frame are satisfied if one sets 

eo'^^Jo^ b= 1,2,3, 



. a 

so that the spatial connection coefficients are given by 

k 

Tbd = {hdbx" - hd^Xb) , 6, c, d = 1, 2, 3, 
with x^ = {lo / 0)5^^ cx'^ . The remaining nonvanishing connection coefficients are 

fbd = Hhbd , l\d = Hh\, b,d= 1,2,3. 
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5.2 Linearised evolution equations 

In this subsection, wc derive the Hnearised system corresponding to the nonhnear equations of Theorem[T] 
for a FRW background with a selfinteracting scalar field. In order to perform the linearisation procedure 
we compute 

and drop all (nonlinear) terms of coupled perturbations. In this way, we obtain the following system 



dt(t> = ■0, 



d^V 



(63) 



\ a / \ 1 h nrn 



2 dV\ „ fdH\ 



V' d^ / \dt J 



ate," = -i/eb" 



a-b, 



1 /UJ 



dtTtd = -^TbdX - Tbp {x''^)/ + H {Srab - hbd&n 
+ BdpeP\-Hrbd- 
As a consequence, the linearised system has the following form 

A^dtii- A.^it,x)djU = B{t,x)u. (64) 
If one defines ii in the same way as in equation (j58p one obtains the linearised matrix B(i,a;) given by 

/B^^^^2 ^2x3 ooooo\ 

-^-BltL 

-4^3x3 213x3 

-it^^Jls (}V' + 2i/)l3x3 ^I3x3 

-^ij^hx3 -^3x3 ^Bi'^^ 

st^fL -H13X3 

I3x3 — i7l3x3 

-^B^xl -^'Bi'xl -fflgxa 

V -jt^Bi'^'l ^Bi'^l -HBi\'l Bi'^l -HI.^J 



where 



r(1) _ / 1 \ (2) _ / U U U \ (3) 

^2x2 — \ i)« _3/f / ' -"2x3 — \ „/, „/, J, ) ' -^3x2 



3^ -giJ 























x^6^ 







x^S^ 


^^6] 







n^: - 






U^' - 






V,v - 






( ^2 


Xj6i 










xjSi 


V 




XjS^ 
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x^S^ 


x^S] 


x'sy 



B 



(8) 
3x3 





xH] 



-x^8] 



x^i] 




B 



(9) 
3x3 



-X^i) 
X^b] 



x^5] 



-x^S^ 



-x^5f 
x^S^ 



B 



(10) 
9x3 



A 


1 


1\ 





























1 


1 


1 





























^1 


1 


V 



B 



(11) 

9x3 



/O 0\ 

1 

1 

1 

1 
1 
1 

\0 0/ 



B 



(12) 
9x3 



fS]x3 S]x^ S]x^\ 

6]x^ S]x^ S]x^ 



6^x1 S]x^ S^x^ 

S]x^ 5^x^ 5]x^ 



SjX^ S^x^ S^x' 

SjX^ SjX^ 5^x^ 

V / 



B 



(13) _ 
9x3 ~ 



/ S^Xj 













S^Xj 












~S]x^ 




Slxj 
















S^x^ 




S^Xj 


-6]x^ 










Sixj 













S^x^ 




\-S^xi 





S{xj 


) 



B 



(14) _ 
9x3 ~ 



/I 














1\ 


1 
















-1 

























1 





B(15) _ 
' ■'-'9x3 ~ 





-1 








1 








1 































1 




1 














1 




-1 








Vo 









^0 








This matrix B is niatic of blocks arounti the diagonal and is close to being of lower triangular form. 
In fact, it could easily be cast into that form by 'removing' the upper matrix B*^^). In the more general 
case where k ^ one has three main blocks. A first block consists of the matrices B^^^, B^^^, B^'^^ B*^''\ 
related to the perturbation variables {(j), ip, xn, X22, X33)- A second block is related to the variables 



iix 



ST- 



bd, dc, 



) and a third one for the perturbation of the spatial connection and the 



remaining spatial frame coefficients, 
results. 



This last block does not play an important role in the stability 



5.3 Results on the nonlinear stability of symmetric hyperbolic systems 

In this subsection we discuss the results of the theory of symmetric hyperbolic systems required to discuss 
the nonlinear stability of perturbations of a reference solution to the Einstein- scalar field system with flat 
spatial sections — i.e. fc = 0. The analysis carried out in the previous sections shows that the Cauchy 
problem for these nonlinear perturbations takes the form 

(AO - eA°{t, u, e)) ^^u - (A^"(i) + eA^{t, u, e)) d,u = (B(t) + eB(t, ii, e)) ii, ^^^^ 
u(a;,0) = Uo(a;), 

where A" is a constant diagonal matrix with positive entries, while A-' and B are matrices whose entries 
are smooth functions of the time coordinate t. The matrices A^, j = 1, 2, 3, are symmetric. Furthermore, 
A°, A-' and B vanish at least linearly with u. If fc = 0, the linearised system ([M)) reduces to 

A°dtu- A^{t)djU = B{t)u. 

Our analysis will require the following definition: 

Definition 1 (Relaxed stability eigenvalue condition). A system of the form given by (|65p is said to 
satisfy the relaxed stability eigenvalue condition if the following conditions hold: 
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(i) There is a constant Sq with < 5q < S{t) sueh that, the eigenvalues \{t,^) of the symbol 

i^jA^'(t)+B(t), i = \/^. 

satisfy 

Re{X{t,i))<-S{t) 
for all n-tuples of integers ^ = (Cij •■•7 ^n) 7^ 0. 

(ii) There is a constant (Sq with < (5o < 5[t) such that, for ^ = 0, the eigenvalues A(t, 0) of^{t), 
satisfy either 

Re {X{t, 0)) < -S{t) or X{t, 0) = 0. 

// the multiplicity of the vanishing eigenvalue is m, then there are m linearly independent eigenvec- 
tors. 

(Hi) The Kernel 0/ B contains the Kernel o/B. 

In order to state the result that allows us to ascertain the nonlinear stability for the solutions of 
the Cauchy problem we introduce some notation. Let H'' {T";R') be the space of all summable 

functions u(-, t) : T" W such that for each multi-index \a\ < k, D°'u{x, t) exists in the weak sense and 
belongs to (T"). The norm in H'' (T" ; R") is defined by 



where 




Given a perturbation u{x,t), we denote by u(^,i) the Fourier coefficient corresponding to the frequency 
vector ^ — recall that it is assumed that the topology of the spatial sections is taken to be that of T'^. In 
particular, u(0,<) denotes the coefficient corresponding to ^ = 0. 

The following theorem shows the exponential decay or convergence to constant values of a solution 
to the nonlinear system (|65p . Its proof follows closely the methods of [HI [H] and will be omitted for the 
sake of brevity. 

Theorem 2. Let P denote the projector into the Kernel o/B(t). Furthermore, let u^*^) = 'Pu{0,t) and 
V = u — u'"' . // the symmetric hyperbolic system of the Cauchy problem (j65p satisfies the relaxed stability 
eigenvalue condition, then the system is stable in the sense that, for sufficiently small e there is eo > e > 
such that 

lim ||v(t)||^n+2(T- = 
and u^°)(t) — >■ const, as time goes to infinity. 

Remark. The use of Fourier methods in Theorem [5] does not allow to disentagle which fields (or 
combination thereof) decay to a constant. However, it is possible to identify the blocks in B(i) responsible 
for this behaviour. 

In the remainder of the article, it will be analysed under which conditions the symmetric hyperbolic 
evolution equations for the perturbations of a reference solution to the Einstein-scalar field system satisfy 
the requirements of the relaxed stability eigenvalue condition (Definition [TJ so that Theorem [5] can be 
used. 
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5.4 Satisfying Assumption (ii) of Definition [T] 



In the case of a flat background (fc = 0, a; = 1), the hnearized matrices A and B are functions of time 
only. Because of the block structure of the matrix, the characteristic polynomial of B is the product 
of the characteristic polynomials of the block matrices around the diagonal. A computation gives the 
following polynomial: 



\ \ ^ J dt ^ 
X ( + bHX"^ - (2il'^ - 6H^ - V") - 




AH' 



(66) 



The first term in the expression ([55)) is the characteristic polynomial of the block related with the per- 
turbation variables {Bhd, e^-', 7"ab)- The second term is related to the variables (fie, ej,°), and the third 
one to {{x^'^)bd, Em)- Finally, the last term — the fifth order polynomial — arises from the block of 
unknowns (0, -0, Xpp). 

In order to obtain conditions from the characteristic polynomial we will make use of results from the 
Lienard- Chipart theorem. The latter gives necessary and sufficient conditions for a polynomial with real 
coefficients to have roots with negative real part — see e.g. [23]. This type of polynomials are called 
Hurwitz polynomials. 

Theorem 3 (Lienard-Chipart). Let 

f{z) = aoz'^ + aiz""^ + 02^""^ H h a„, oq > 0, 

be a polynomial with real coefficients. Then the following statements are equivalent: 
(i) the polynomial is a Hurwitz polynomial; 

(ii) The coefficients of f are positive and 62 > 0, di> 0, ...,(5„, n even; 
(Hi) The coefficients of f are positive and Si > 0, 63, > 0, n odd, 
where the Hurwitz determinants are defined by 





/ fli 


03 


05 • 


• a2i-i\ 




1 


a2 




O21-2 


So = 1, Si = det 





ai 


0.3 ■ 


021-3 




U 





• 


oi J 



It is easy to see that the first term in (|66p requires 

H{t) > 



(67) 



which is simply the condition for an (ever) expanding background. From the second factor in the char- 
acteristic polynomial (|66p. one obtains the following conditions on V': 



V H 



2hZ<-^-2H^ 
iP dt 



Notice that by construction we have that i/j < — see equation ([5^ . The first condition in (|68p implies 
that the first derivative of the scalar field potential is positive. The second condition is related to the 
evolution equation for the Hubble function (pTj) . Moreover, by virtue of condition ([50)) the ratio V'/i/j 
tends to a constant which is bigger than the asymptotic value of H/2. 



The fifth order polynomial in (j66|) can be rewritten as 

A2 (x 



5HX^ 



2V'^ - 6H^ - V" ) A + i^2HV'' 
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from where one sees that the matrix B (t) has a double vanishing eigenvalue. A direct application of the 
Lienard-Chipart Theorem to the remaining third order polynomial renders the following two conditions: 

The first condition relates the second derivative of the potential to its first derivative, and implies that 
V" has to be positive and bounded. The second condition is sharper and also relates to the evolution of 
the Hubble function (pT|) . Taking the limit ?/! — > 0, one recovers the necessary asymptotic conditions for 
stability found in [11) . 



5.4.1 Nonflat backgrounds 

For the sake of conciseness, the case of a reference solution with nonfiat spatial sections will not be treated 
here. Nevertheless, we point out some features of the required analysis. For this class of solutions the 
characteristic polynomial of the matrix B is given by 

/(A) X (^X^ + 5HX^ - (2ij^ - 6H^ - V") A''' + (2HV" + A^^ x (A + (70) 

where /(A) is a polynomial of degree 15 in A. The coefficients of this polynomial are, unfortunately, 
large expressions without an obvious structure. Accordingly, an analysis of the roots by means of the 
Lienard-Chipart Theorem, although in principle possible, becomes extremely cumbersome. The study 
of perturbations of cosmological models with nonflat spatial sections also requires to take into account 
extra terms coming from partial differentiation when constructing the Sobolev norms. This is because 
the matrices A-' and B are now functions of the coordinates . 

Finally, it is mentioned that a discussion of the perturbations of reference solutions with nonflat spatial 
sections requires suitable adaptations of the relaxed stability eigenvalue condition (Definition [1} and of 
Theorem H 



5.5 Satisfying Assumption (i) of Definition [T] 

The purpose of the present section is to show that Assumption (i) of the stability eigenvalue condition 
follows from Assumption (ii) . This turns out to be a consequence of a theorem due to Ortiz in [21] — see 
also [inj. 

Theorem 4. Let S be a nonsingular matrix which takes the s x s matrix B into block diagonal form. If 
for all unit frequency vectors ^, the linear map 

iS-iA°S + iA(i,|) : Ker{B{t))^C' 

with 

A(i,0 = -iS-i^,A^S, 

is injective, then there are positive constants o-nd fxg — 6q/c, such that all eigenvalues A(t,^) of the 

symbol i^jAP (t) + B(t) satisfy 



i?eA(t,0 < 



-fiom^ or - So, if III <eo 
--50, if \i\ > ^0 



Remark. This theorem states that if the determinant of the linear map is different from zero, then the 
eigenvalues of the linearised symbol have negative real part for nonzero Fourier frequency, decaying no 
faster than ||p when ||| 0. 

In order to apply this theorem, we must first put B into block diagonal form. To this end we need to 
find a transformation, S, such that 
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To achieve this, it is sufficient to consider the block matrix corresponding to the degree 5 polynomial 
which contains the two zero eigenvalues. Let us define the 5x5 matrix M by 



M 



^r(I) r(2) ^ 
■'-'2x2 ■'-'2x3 



B 



(3) 
3x2 



B 



(4) 
3x2 / 





1 










-V" 


-2,H 


-V' 








-l^ 










-It 










-¥ 









Once M has put into block diagonal form, it turns out that the condition ([TT]) holds. In order to find 
a suitable transformation, S, we first make a change of basis taking M into block triangular form. For 
this, we compute the eigenvectors corresponding to the two vanishing eigenvalues. From the condition 



Mv = 0, 

it follows that the eigenvectors can be taken to be 



/ V" 1 



Vl 




1 
1 

V 1 J 



V2 




1 
1 



With the help of these two eigenvectors we can easily transform M into block triangular form by com- 
pleting the basis with three further linearly independent eigenvectors. For instance, one can consider 



V3 = 



1 





V4 = 





1 



VV 



V5 = 






1 

VV 



Using these basis vectors as the columns of a transformation matrix one arrives to 





3i> 

V" 


&H 
V 








o\ 










1 










1 


1 





1 







1 


1 








1 


V 


1 


1 





1 


1/ 



with inverse given by 



v"v' 



3 2HV"+i)V' 
1 V"V'_ _ 
3 2HV"+iliV' 










1 






2gV" _ 
2HV"+i<V' 

2HV"+i<V' 




-1 



2HV" _ 
2HV"+i>V' 

2HV"+i<V' 



-1 




2HV" 
' 2HV"+il>V' 

' 2HV"+il)V' 



1 
1 



As a result of these computations one finds that 



P^MP = 



L3x3 
02x3 



H3X2 
02x2 



with 



-■3x3 



/ -2H 







2HV' 



V" 4_H"V'.+V' 
3 2HV"+',pV 

V' 2V)^-V"^ 
3 2HV"+i!V' 

-3H 



H 



■3x2 



3 2HV"+ipV' 

4 -H"V'V'_ _ 
3 2HV"+i!V' 



8 H^V'[ \ 
3 2HV"+'ipV' 

' 3 2HV"+ipV' 

-2i, J 
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It can be readily verified that the matrix L is negative definite as long as the conditions of the 
previous section hold. Now, in order to proceed further we need another transformation which converts 
the block upper triangular form into block diagonal. Choosing 



R: 



I3x3 
02x3 



Q3x2 
T-2X2 



R 



one finds that 



R-l (P iMPl R 



I3x3 
02x3 

LQ + H 




-Q3x2 
I2x2 



Thus, in order to obtain the required block-diagonal form we just need to solve 



A computation shows that 



LQ 



HV" 



H. 



HV" 



Q = 



V 



3 2HV" + ipV' 

2 4>V' 

3 2HV"+'ipV' 





3 2HV"+i>V' 

2 _ 4>V' 

3 2HV"+ipV' 





so that the whole transformation matrix S is given by 



where 



PR 



(PR)5X5 
030 x5 



05x30 
130x30 



/- 


3'4> 


6H 





V" 


V 










1 




1 


1 







1 


1 





V 


1 


1 







1 

1^ 

3 



_ 2 
1^ 

! 

3 ' 



with inverse given by 



(PR) 



1 v'v" 





1 


2 HV" 


2 HV" 


3 2HV"+i>V' 
1 V'V" 


3 2HV"+4'V' 
1 JpV' 


3 2HV"+4:V' 
1 i/>V' 


3 2HV"+i>V' 




3 2HV"+4;V' 




3 2HV"+4:V' 





















-1 






HV" 



3 2HV"+i>V' 

1 rpV' 

3 2HV"+i/jV' 



1 
1 



Summarising: we have constructed the matrix S which transforms the linearised matrix B into block 
diagonal form. In the following we use this matrix to show that the linear map arising in the formulation 
of Theorem [His injective. 



Applying the transformation S, to the linearised matrices of equation (|59|) — see also equation ((63 
we get that 

/ 




where 



Ad) 



2w 



HV" 



3a 2HV"+ipV 
3a 2HV"+4>V' 






7^1 

76 



J3) 
^-3x5 









HV" 






(A(^))L3 






,(1) 



-^5x3 
a(2) 
-^3x3 





3a 2HV"+il:V' 
UJ ^V' 



3a 2HV"+-4:V 




76 
76 



2ui 



HV" 







(A(^))Jx3 




3a 2HV"+-iliV 



76\ 



3a 2HV"+4>V 





76 



a?3 

UJ f 

a?3 



A(2 






o\ 


























-^3x3 






















0^ 










UJ c 
a?l 







^?3 












UJ f 

a« 


^6 
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It is also noticed that the diagonahty of A" imphes that it is invariant by S. It fohows that the hnear 
map iA" + iA(i, ^) has determinant 



where = + ^| + Therefore, for aU ^, a and uj the determinant is nonzero, and by Theorem |4] 
one finds that Assumption (i) of Definition [1] holds. 



5.6 Verifying Assumption ('mj of Definition [T] 

In order to verify Assumption (in) of the relaxed stability eigenvalue condition, we just need to consider the 
block matrix of the nonlinear matrix, B, associated to the perturbation variables (0, ip, xn? X22, Xaa)- 
A simple inspection of the linearisation process allow us to verify that this block matrix has zeros in all 
its entries. Therefore Assumption (Hi) follows directly. 



6 The main result 

We summarise the discussion of previous sections in the following theorem: 

Theorem 5. Let (f> be a smooth real scalar field with a selfinteracting potential V((/)) in a flat expanding 
Friedmann- Robertson- Walker spacetime subject to the nonlinear evolution equations 

d0 diA ; dV 

— — = %h and — — = —?>Hib -. 

At ^ At ^ 

Suppose there exists a constant Hq such that < Hq < H{t), for all t > 0, and a potential satisfying 
conditions 

Then the FRW-nonlinear scalar field solution is stable in the sense that, given any initial data for small 
nonlinear perturbations Uq whose norm ||'u|j//.">(T3) is finite, there exists eo > such that for all < e < sq 
a solution to the nonlinear perturbations exists for all times and decays exponentially to zero or converges 
to constant values. 



Remark. Our main result can be understood as a cosmic no-hair theorem for a wide class of scalar 
field cosmologies. We note that our results can be useful when considering numerical simulations of these 
models as it gives conditions for exponential decay during the evolutions. It could, therefore, be used to 
test the consistency of the numerical steps. 



6.1 Power law solutions with exponential potential 

As a direct application of Theorem [5] we discuss the so-called power-law solutions to an exponential 
potential and show that this class of solutions is covered by our analysis. 

Consider the reference solution for a scalar field (p with exponential potential 

V((/.) = Ae^^ 
with A, A positive constants, given by the metric 

ds2 = -At^ +t^P5i,dx'dx^ 
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with 




with p > 1/3 and X^p = 2. A direct computation shows that the conditions in Theorem [S] are satisfied if 

p > « 0.454. 

Thus, for this range of values one can conclude that sufficiently small nonlinear perturbations to power-law 
solutions decay exponentially to zero or converge to constant values. 

Remark 1. It is well known that for p > 1 (i.e. A < a/2), the exact solution described in the previous lines 
has accelerated expansion — the so-called power law inflation. Our main result. Theorem [S] shows that 
one has decay of the nonlinear perturbations even in the absence of accelerated expansion if p g ( ^^^^ ; !]■ 
This is, to the best of our knowledge, a new result — the possibility of having decay in the marginal case 
p = 1 has been discussed in [T2j [T3| . 

Remark 2. Simple numerical simulations of the evolution of linearised perturbations supports the results 
of our analysis. A detailed discussion of these is beyond the purposes of the present work. 
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